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Abstract 

We consider a multidimensional extremal problem formulated in 
terms of tropical mathematics. The problem is to minimize a nonlinear 
objective function, which is defined on a finite-dimensional semimodule 
over an idempotent semifield, subject to linear inequality constraints. 
An efficient solution approach is developed which reduces the problem 
to that of solving a linear inequality with an extended set of unknown 
variables. We use the approach to obtain a complete solution to the 
problem in a closed form under quite general assumptions. To illus- 
trate the obtained results, a two-dimensional problem is examined and 
its numerical solution is given. 
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1 Introduction 

Models and methods of tropical (idempotent) mathematics [U El [31 HI El [7] 
find expanding applications in many fields including engineering, computer 
science and operations research. There are a variety of real-life problems that 
can be represented and solved as tropical extremal problems. As related 
examples one can take idempotent algebra based solutions to scheduling 
[TJ U M U ED] and location [2 QU [El Q31 Q31 US E] problems. 

The tropical extremal problems are usually formulated as to minimize 
functionals, which are defined on finite-dimensional semimodules over idem- 
potent semifields, subject to linear constraints in the form of equalities or 
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inequalities. The problems under study may have linear objective functions 
as in pH El [7] or nonlinear function as in [T^ [T^ [T^ [2QI [131 [Ml [13 Ell [J6l [22] . 
In some cases both the objective function and the constraints in the prob- 
lems appear to be nonlinear (see, e.g. [T2"l 19]). 

Many existing solution techniques [T2J El [3 [21] offer an iterative al- 
gorithm that gives particular solutions if any, and indicates that there is no 
solution otherwise. Other approaches [III [JTl [121 EDI [13 [111 [151 [13 122] allow 
one to get direct solutions in a closed form. However, in most cases only 
partial solutions are given leaving the problems without complete solution. 

In this paper, we consider a multidimensional problem with a non- 
linear objective function under linear inequality constraints. The prob- 
lem is actually a constrained extension of the unconstrained problems in 
[T71 \19\ [20] \13\ [HI [IS]- A solution approach is developed which reduces the 
problem to solving a linear inequality with an extended set of unknown vari- 
ables. We use the approach to obtain a complete solution to the problem in 
a closed form under quite general assumptions. 

First we outline basic algebraic facts and present preliminary results 
that underlie subsequent developments. Together with an overview of known 
results including the solution of linear equations, a new binomial identity for 
matrix traces is presented. Furthermore, we formulate the extremal problem 
of interest and then derive a complete solution. To illustrate the obtained 
results, a two-dimensional problem is examined and its numerical solution 
is given. 

2 Preliminary Results 

The aim of this section is to provide a basis for further analysis and solution 
of tropical extremal problems. We briefly outline key algebraic facts and 
present preliminary results from |181 12U] . For further details and considera- 
tions, one can refer top^EllIlElElE]. 

2.1 Idempotent Semifield 

Let X be a set that is endowed with addition © and multiplication ® and 
contains the zero element and the identity 1. We suppose that both 
operations are associative and commutative, and multiplication is distribu- 
tive over addition. Moreover, addition is idempotent, which means that 
x (& x = x for all x € X , and multiplication is invertible; that is there exist 
an inverse x _1 for any x 6 X + , where X + = X\{0}. Due to the above prop- 
erties, (X, 0, 1, ©, <g>) is usually referred to as the idempotent commutative 
semifield. 

The power notation with integer exponents is routinely defined in the 
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semifield. For any x G X + and an integer p > 1 , we have x° = t, (D p = , 
x p = x p_1 <S)X = x®x p ~ l , x~ p = (x~ 1 ) p . 

Moreover, we assume that the rational powers are also defined which 
makes the semifield radicable. 

In what follows, the multiplication sign (g> is omitted as is common in 
conventional algebra. The power notation is used only in the above sense. 

Idempotent addition induces a partial order on X so that x < y if and 
only if x ® y = y . We assume that the partial order can always be extended 
to a total order and so consider the semifield as linearly ordered. From here 
on, the relation symbols and the operator min are thought of in terms of 
this linear order. 

An example of a commutative idempotent semifield that is radicable and 
linearly ordered is the real semifield R m ax,+ = (RU {— oo}, — oo, 0, max, +} . 



2.2 Idempotent Semimodule 

Consider the Cartesian product X n with column vectors as its elements. For 
any vectors x = (xj) and y = (yi) from X n , and a scalar c G X, vector 
addition and scalar multiplication are routinely defined component-wise 



{x®y}i = Xi ®y u {cx}i 



CXi • 



A geometric interpretation of the operations for the semimodule R 
is given in the Cartesian coordinates in Fig. [TJ 



2 

max, 




CX 2 



yi xi 




xi cxi 



Figure 1: Vector addition (left) and scalar multiplication (right) in K max + . 

Equipped with these operations, the set X n forms a semimodule over 
the idempotent semifield X . 

A vector with all zero elements is called the zero vector. A vector is 
regular if it has no zero elements. The set of all regular vectors in X n is 
denoted by X!L . 
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For any nonzero column vector x = (xj) , we define a row vector x = 
(x~) with the elements x~ = x~ x if x% > 0, and xj = otherwise. 

For any regular vectors x and y , the component- wise inequality x < y 
implies x > y . 

A vector y is linearly dependent on vectors x\, ... ,x m if y = c\X\ © 
• • • c m x m for some scalars c±, . . . ,c m . Specifically, vectors y and x are 
collinear if there exists a scalar c such that y = ex. 

Given vectors x±, . . . ,x m , the set of linear combinations c\X\ • • • © 
c m x m for all ci, . . . , c m is a linear span of the vectors, which forms a sub- 
semimodule. A graphical example of the linear span of two vectors is shown 
in Fig. [21 
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Figure 2: Linear span of two vectors in 



2.3 Matrix Algebra 

Let X nxn be the set of square matrices of order n with entries in X . For any 
matrices A = (ciij) and B = (bij) from X nxn , and a scalar c € X, matrix 
addition and multiplication together with scalar multiplication follow the 
conventional rules 

n 

{A © B}ij = dij © bij, {AB}ij = ^^a-ikbkji {cA}ij = caij. 

fc=i 

A matrix with all zero entries is called the zero matrix and denoted by 
(D . A matrix with all off-diagonal entries equal to is a diagonal matrix. If 
all entries of a matrix above or below the diagonal are zero, the matrix is 
triangular. A diagonal matrix with 1 on the diagonal is the identity matrix 
denoted by /. 

A matrix is reducible if it can be put into a block-triangular form by 
simultaneous permutations of rows and columns. Otherwise, the matrix is 
irreducible. 
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The matrix power is introduced in the regular way. For any matrix A 
and integer p > 1 , we have 

A° = I, A p = A p ' 1 A = AA p ' 1 . 

2.3.1 Matrix Trace 

For any matrix A, the trace of A is routinely given by 

tr A = an © ■ ■ ■ © a nn . 

It is not difficult to verify that for any matrices A and B , and a scalar 
c it holds 

ti{A®B) = tr A®txB, tr(AB) = tr(BA), 
tr(cA) = ctr(A). 

We apply the above properties to derive a binomial identity for traces 
to be used below. 

Lemma 1. For any matrices A and B , and an integer m > 1 it holds that 

m 

tv{A®B) m = tvB m ® tr (AB h ■■■ AB ik ). (1) 

k=l iiH \-i^=m— k 

Proof. First we note that the power (A © B) m can naturally be expanded 
as the sum of 2 m products formed by k matrices A and m — k matrices B 
taken in all possible orders for all k = 0, 1, . . . m. 

Furthermore, we take trace of the sum and consider those terms that 
have at least one matrix A as a factor. Since two matrices commute under 
the trace operator, all these terms can be rearranged so that each term has 
A as the first factor. It remains to see that any product with k > 1 matrices 
A is then represented as (AB n ) ■ ■ ■ (AB lk ) , where are nonnegative 

integers such that i\ + ■ ■ ■ + = m — k . □ 

2.3.2 Spectral Radius 

As usual, a scalar A is an eigenvalue of a matrix A, if there exists a nonzero 
vector x such that 

Ax = \x. 

Every irreducible matrix has only one eigenvalue, whereas reducible ma- 
trices may have several eigenvalues. The maximal eigenvalue (in the sense 
of the linear order on X) is called the spectral radius of A. It is directly 
calculated as 

n 

A= 0tr 1 /" t (^ m )- 

m=l 
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The spectral radius A of a matrix A offers a useful extremal property, 
which holds that 

min x~ Ax = A, 
where the minimum is over all regular vectors x . 

2.4 Linear Inequalities 

Let A £ x nxn be a given matrix, x € X n be the unknown vector. Consider 
the problem of finding regular solutions x to the inequality 

Ax < x. (2) 

For any matrix A , the solution of the inequality involves a function that 
is given by 

Tr(,4) = tiAe---etrA n , 
and a star operator that is defined as 

A* = I ® A® ■ ■ ■ ® A™" 1 . 

In the case of irreducible matrices, the solution is given by the following 
result. 

Theorem 2. Let x be the general regular solution of inequality ([2]) with an 
irreducible matrix A. 

Then the following statements hold: 

1. If Ti{A) < 1, then x = A*u for all u € X™ . 

2. If Tr(A) > 1, then there is no regular solution. 

Fig. offers an example of solution of a linear inequality with a matrix 

2 

max,+ • 



A = (ai,a 2 ) in R 2 



3 A Constrained Extremal Problem 

Given matrices A, B € X nxn , the problem is to find all regular solutions 
x € X™ so as to provide 

min x Ax, 
Bx < x. 

Below we give a general solution to the problem and then illustrate the 
solution with numerical examples in the framework of the semifield (R max + • 
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Figure 3: Solution to a linear inequality in 



3.1 The Main Result 

A complete direct solution to the problem under quite general conditions is 
given as follows. 

Theorem 3. Suppose that at least one of matrices A and B is irreducible, 
A > is a spectral radius of the matrix A, and Tr(B) < t . 
Then the minimum in Q is equal to 

n 

= tr 1 /* (AB* 1 ■■■AB^) (4) 

k=l 0<iiH Mfc<n— k 

and it is attained if and only if 

x = (6- l A@B)*u (5) 

for all regular vectors u . 

Proof. First we introduce an additional variable and then reduce problem 
([3]) to an inequality. Existence conditions for solutions of the inequality is 
used to evaluate the variable, whereas the solution of the inequality is taken 
as the solution of the problem. 

Suppose 9 is the minimum of the objective function in problem ([3]) and 
note that 9 > A > . 

The set of regular vectors x that yield the minimum is determined by 
the system 

x Ax = 9, 
Bx < x. 

Let us examine the first equality x~ Ax = 9. Since x Ax > 9 for all 
x € X!f , the equality can be replaced by the inequality x Ax < 9. 
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Furthermore, we multiply the inequality by O x from the left. Since 
x x = 1 and xx~ > I for all regular x, we have 0~ l Ax < Q~ x xx~ Ax < x 
and so get the inequality _1 Ax < x. Considering that left multiplication 
of the last inequality by 0x~ gives the first one, both inequalities prove to 
be equivalent. 

Now the solution set of the problem is given by the system of inequalities 

0~ x Ax < x, 
Bx < x. 

The above system is equivalent to one inequality 

(0~ X A © B)x < x, (6) 

which takes the form ([2]) with an irreducible matrix. 

It follows from Theorem [2] that inequality (|6|) has regular solutions if and 
only if 

Tr(0 -1 jl © B) < 1. (7) 
Consider the function Ti(0~ l A © B) and represent it in the form 

n 

Tr(Q- l A © B) = triO' 1 A © B) m . 

m=l 

By applying binomial identity (JTJ to the right-hand side, we have 

n m 

TriO^AeB) = TrSffi 0~ k ix{AB h ■ ■ ■ AB**). 

m=l k=l ixH hifc=m— k 

Furthermore, rearrangement of terms in the last sum gives the following 
expression 

n 

Tr{0~ l A®B) = Tr5ffi0 0~ k tr(AB n • • • AB ik ). 

k=l 0<nH Uk<n-k 

Considering that Tr B < 1 by the conditions of the theorem, inequality 
d?]) reduces to inequalities 

0~ k tr(AB h ■■■AB lk ) < 1, 

0<«iH \-i k <n— k 

which must be valid for all k = 1, . . . , n . 

By solving the inequalities with respect to , we get inequalities 

0> ti 1/k (AB h ■■■AB ik ), 

0<uH \-ifc<n—k 
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or equivalently, one inequality 

n 

^-0 © tr 1 / fc (AB il ---AB i *). 

k=l 0<«iH \-i k <n-k 

In order for 9 to be the minimum in problem ([3]), the last inequality 
must be satisfied as an equality, which gives (jU). 

It remains to apply Theorem [2] to inequality Q so as to arrive at the 
solution in the form ([5|). □ 

As a particular case of the result, a solution to problem [3] without linear 
constraints can be derived. 

Corollary 4. Let, in addition to the assumptions of Theorem 0, A be an 
irreducible matrix and B = . 

Then the minimum in ([3D is equal to A and it is attained if and only if 

x = (\~ l A)*u 

for all regular vectors u . 

Proof. It is sufficient to verify that we have 8 = A. For this purpose we 
rewrite dH in the form 

n n—1 

= 0tr 1 / fc (A*)©0 tr 1 1 k {AB h ■■■AB^). 

k=l k=l 1 <ii H hifc<«— k 

With B = we arrive at the desired result. □ 
Note that this result is consistent with that in |22j . 

3.2 Illustrative Examples 

Consider problem ([3]) in R^ ax + with matrices 

First we apply Corollary [3] to solve the problem without constraints. We 
calculate the matrix 



(2 



-3 
-5 -4 



and then find 

A = tr^©tr 1/2 (^ 2 ) = = 1. 
Since A _1 ^4 = A, we have 

(X^A)* = A* =I®A = ( _° 5 ~ 3 Q 
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The solution to the problem takes the form 

-3 



x 



-5 



Note that in a similar way, one can obtain the solution to the inequality 
Bx < x in the form 



x 



A graphical illustration of the solutions on a plane with the Cartesian 
coordinates is given in Fig. 21 




Figure 4: Solutions to an unconstrained problem (left) and inequality con- 
straints (right) in R^ax + • 

Now we use Theorem [3] to handle the constrained problem. We get the 
matrix 

and then calculate 

9 = tr A © tr 1/2 (A 2 ) © tr(AB) = 2. 
Furthermore, we have 

and then arrive at 

' -5 



{e~ 1 A®By 



5 o 



Since the columns in the obtained matrix are collinear, we take only one 
of them, say the first. The solution set is then given by 





X = \ ) U = UXq 
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for any nonzero number u G ^max,+ • 

Fig. [5] combines the solutions to the unconstrained problem and to the 
constraints with that of the constrained problem. The solution of the last 
problem is depicted with a double-thick line that is drawn through the end 
point of the vector x$ and coincides with a border of the feasible area defined 
by the constraints. 




Figure 5: Solution (double-thick line) to the problem under inequality con- 
straints in R?.„ , . 
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